Extinction of quasiparticle scattering interference in cuprate superconductors 
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The quasiparticle scattering interference phenomenon characterized by the peaks in the local 
density of states is studied within the kinetic energy driven superconducting mechanism in the 
presence of a single impurity. By calculation of the Fourier transformed ratio of the local density 
of states at opposite energy, it is shown that the quasiparticle scattering interference phenomenon 
can be described qualitatively by a single impurity in the kinetic energy driven homogeneous d- 
wave superconducting state. The amplitude of the peak increases with increasing energy at the low 
energy, and reaches a maximum at the intermediate energy, then diminishes to zero at the high 
energy. The theory also predicts that with increasing doping, the position of the peak along the 
nodal direction moves towards to the center of the Brillouin zone, while the position of the peak 
along the antinodal direction is shifted to large momentum region. 
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In the study of the physical properties of cuprate super- 
conductors, one of powerful techniques is scanning tun- 
neling microscopy (STM)^, since it is the only method 
to probe the real-space inhomogeneous electronic struc- 
ture of cuprate superconductors. Application of STM 
to cuprate superconductors would allow one to explore 
some important issues^^^, including (1) the physical pro- 
cesses dominating quasiparticle scattering, (2) the degree 
to which the quasiparticles are well defined and coherent, 
and (3) the relation between the commensurate and in- 
commensurate magnetic signatures and the quasiparticle 
scattering processes. In particular, STM has been used to 
infer the momentum-space information of quasiparticles 
from the Fourier transform of the position (r) and energy 
(a;)-dependent local density of states (LDOS), p(r,w), 
then both real-space and momentum-space modulations 
for LDOS are explored simultaneously^. The typical fea- 
ture observed by the Fourier transformed LDOS is dom- 
inated by the peaks at well-defined wavevectors obey- 
ing the octet modeP"^, since the quasiparticle disper- 
sion in cuprate superconductors in the superconducting 
(SC) state has closed constant-energy contours surround- 
ing the d-wave nodes. Theoretically, there is a general 
consensus that the peak emerges due to the quasiparticle 
scattering interference (QSI)^"— liS. This QSI is mani- 
fested itself as spatial modulation in the Fourier trans- 
formed LDOS /o(q, w) whose wavevectors change with en- 
ergy. Several attempts have been made to make this ar- 
gument more precise. In particular, some calculations 
based on the phenomenological d-wave Bardeen-Cooper- 
Schrieffer (BCS) formalism have been performed by con- 
sidering the effect of impurity scatteringS^--, where it has 
been shown that a single or few impurities in a homoge- 
neous d-wave SC state leads to a result that is in quali- 
tative agreement with the STM experimental data-"—. 

Recently, the improvements in the resolution of STM 
experimentsi^rJi allowed one to resolve additional fea- 
tures in QSI. Among these new achievements is the mea- 
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surement of the ratio of differential conductances at op- 
posite bias, 

Z{r,V)=g{r,+V)/g{r,~V), (1) 

where V is the bias voltage and g{r, V) the differential 
conductance. The advantage of this procedure is that 
it cancels the severe systematic errors in the measure- 
ment of g{r, +V) due to tip elevation errorsiSrJl, yet re- 
tains all the quasiparticle information in the differential 
conductance. In particular, an remarkable phenomenon 
observed from these new STM experiments is that QSI 
always disappears at a scale energy (high energy) that is 
indistinguishable from the energy at which electronic ho- 
mogeneity is lost. In other words, at energies below the 
high energy, the STM response is similar to the dispers- 
ing Bogoliubov quasiparticle, however, at energies above 
the high energy, the response becomes highly spatially in- 
homogeneous. In corresponding to the constant-energy 
contour of the high energy, an extinction line exists^'* in 
momentum space [near the diagonal line between (±7r, 0) 
and (0, ±7r)], beyond which most of the dispersing peaks 
disappear, to be replaced by a reduced set of roughly non- 
dispersive peaks. These new STM experimental results 
lead naturally to speculation about the appearance of 
the QSI extinction in cuprate superconductors in the SC 
state near the boundary of the Brillouin zone (antinodal 
regime). Within the phenomenological BCS formalism, 
the QSI phenomenon in a d-wave SC state with coexist- 
ing short-range antiferromagnetic (AF) order has been 
studie d"'"- ' recently by considering the scattering aris- 
ing from a single point-like impurity, and the result shows 
that dispersing peaks are then extinguished at the high 
energy!^. To the best of our knowledge, the dispers- 
ing peaks below the high energy in cuprate superconduc- 
tors and its extinction at the high energy has not been 
treated starting from a microscopic SC theory, and no 
explicit calculations of the doping dependence of the dis- 
persing peaks has been made so far. In this paper we 
show very clearly if the effect of a single point-like impu- 
rity scattering potential is considered within the kinetic 
energy driven SC mechanismii, the calculated ratio of 
LDOS at opposite energy in momentum-space can re- 
produce some main features observed experimentally on 
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cuprate superconductors^i^ii^, including the appearance 
of the dispersing peaks below the high energy and its 
extinction at the high energy. With increasing doping 
concentration, the position of the peak along the nodal 
direction moves towards to the center of the Brillouin 
zone, while the position of the peak along the antinodal 
direction is shifted to large momentum region. 

In cuprate superconductors, the single common feature 
is the presence of the two-dimensional Cu02 plane, and 
then it is believed that the unconventional physical prop- 
erties of cuprate superconductors is closely related to the 
doped Cu02 planes^*. It has been argued that the t-J 
model on a square lattice. 
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acting on the Hilbert subspace with no double occupied 
site, i.e., Cj^Cia < 1, captures the essential physics of 
the doped Cu02 plane^^, where r) = ztx, ±y, t = zLxzty, 
C}^ {Cia) is the electron creation (annihilation) operator. 
Si = {Sf , Sf , Sf) are spin operators, and /i is the chem- 
ical potential. To deal with the constraint of no dou- 
ble occupancy in analytical calculations, the charge-spin 
separation (CSS) fermion-spin theory2£ has been devel- 
oped, where the constrained electron operators are de- 
coupled as Cif — hl^S^ and Cii — ^li^t ^ ^ifh the 
spinful fermion operator hia = e^**'"/ii describes the 
charge degree of freedom together with some effects of 
the spin configuration rearrangements due to the pres- 
ence of the doped hole itself (charge carrier), while the 
spin operator Si describes the spin degree of freedom 
(spin), then the electron local constraint for single oc- 
cupancy, Y.a'^lad'y = Sfhi^h\^Si + S~hi^h\^Sf = 
h,h\{SfS~ + S^S^) = 1 - h\h^ < 1, is satisfied in ana- 
lytical calculations. In particular, it has been shown^^ 
that under the decoupling scheme, this CSS fermion- 
spin representation is a natural representation of the con- 
strained electron defined in the Hilbert subspace without 
double electron occupancy. Furthermore, these charge 
carrier and spin are gauge invariant22i2i, and in this 
sense they are real and can be interpreted as physical 
excitations^^. This is much different from the usual slave- 
boson approach22i2i, where the electron operator is de- 
composed as the holon and spinon, however, the local 
constraint for the single occupancy is explicitly replaced 
by a global constraint in the mean-field level. Due to the 
constraint, these holon and spinon are also coupled by a 
strong gauge fieldS^i^, allowed by this slave-boson repre- 
sentation, and therefore these holon and spinon are not 
gauge invariant. From this point of view, our treatment 
of constraint for the physical electron may be better than 
the usual slave-boson approach. In this CSS fermion-spin 
representation, the t-J Hamiltonian ^ can be expressed 
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with Jeff = (1 — (5)^ J, and 5 = {h\^hia) = {h\hi) is the 
hole doping concentration. As a consequence, the kinetic 
energy term in the t-J model has been transferred as an 
interaction between charge carriers and spins, which re- 
flects that even the kinetic energy term in the t-J Hamil- 
tonian has a strong Coulombic contribution due to the 
restriction of no double electron occupancy of a given 
site. 

For the understanding of the physical properties of 
cuprate superconductors in the SC state, we have de- 
veloped a kinetic energy driven SC mechanism~, where 
the interaction between charge carriers and spins arising 
directly from the kinetic energy term in the t-J model ([3]) 
induces a d-wave charge carrier pairing state by exchang- 
ing spin excitations in the higher power of the hole doping 
concentration, then the electron Cooper pairs originat- 
ing from the charge carrier pairing state are due to the 
charge-spin recombination, and their condensation re- 
veals the SC ground-state. In particular, this d-wave SC 
state is controlled by both the SC gap function and the 
quasiparticle coherent weight, which leads to a fact that 
the maximal SC transition temperature occurs around 
the optimal doping, and then decreases in both under- 
doped and overdoped regimes. Furthermore, it has been 
shown that this SC state is a conventional BCS-like with 
the d-wave synimetry^^, so that the basic BCS formal- 
ism with a d-wave SC gap function is still valid in quan- 
titatively reproducing some main low energy features of 
the SC coherence of quasiparticles, although the pairing 
mechanism is driven by the kinetic energy by exchanging 
spin excitations. Within this kinetic energy driven su- 
perconductivity, we have discussed the low energy elec- 
tronic structure^i^^ of cuprate superconductors, the dy- 
namical spin response^ and the quasiparticle trans- 
port in the SC stated, and qualitatively reproduced some 
main features of ARPES experiments"'^-, inelastic neutron 
scattering experiments^*'^^, and microwave conductivity 
measurements'^'^ on cuprate superconductors. The typi- 
cal feature of this kinetic energy driven SC mechanism 
is that the pairing comes out from the kinetic energy 
by exchanging spin excitations and is not driven by the 
magnetic superexchange interaction as in the resonant 
valence bond type theories In particular, a possible 
SC theory has been developed for description of super- 
conductivity in doped hexaboridea^, where the physical 
mechanism favoring such a reorientation is the enhanced 
coherence (and hence lower kinetic energy) of the doped 
electrons in a ferromagnetic background relative to the 
paramagnet. Following the previous discussion a^^'^^ , the 
charge carrier diagonal and off-diagonal Green's func- 
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tions of the t-J model ([3]) can be obtained as, 
/ Tj2 y2 \ 

52i(k,w) = t; , , , p h (5) 

where the charge carrier quasiparticle spectrum E^k ~ 

\J I A;i2(k) 1 2 with the renormahzed d-wave charge 

carrier pair gap function Ahzi}^) = ^hz[coskx—cos.ky\/2, 
and the charge carrier quasiparticle coherence factors 
UL = (1 + 6c/i?hk)/2 and V^^ = (1 - 6c/i?hk)/2, 
while the charge carrier quasiparticle coherent weight 
ZhF and other notations are defined as same as in Ref. 
|25| . and have been determined by the self-consistent 
calculationi2i2^. For the convenience of the following dis- 
cussions, the full charge carrier Green functions (|4]) and 
([5]) can also be expressed in the Nambu representation 



as, 
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where tq is the unit matrix, ti and T3 are Pauli matrices. 

In the CSS fermion-spin representatio n^°i^^ , the elec- 
tron Green's function in the Nambu representation, 

r(V,,\-( G'ii(k,w), Gi2(k,a;)\ 

is a convolution of the spin Green's function and charge 
carrier Green's function ([5]), and its diagonal and off- 
diagonal components Gii(i—j, = ((Cio-(i); Cj^(t'))) 
and G2i(i -j,t- t') = {{Cl^{t)] C]^{t'))) have been ob- 
tained as^^. 
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respectively, where the electron quasiparticle coherent 
weight Zp — ZhF/'2, and the spin excitation spectrum 
Wp and Bp have been given in Ref. [2^. 

In the presence of a single point-like impurity scatter- 
ing potential, 

V = VoSir)T3, (9) 

the unperturbed electron Green's function in Eq. ([7]) 
is dressed by this impurity scattering, where the T ma- 
trix exactly accounts for multiple scattering off that 
impurity^^. Since translational invariance is broken by 
the impurity, the dressed electron Green's function in 
the Nambu representation in real-space depends on two 
positions r and r', 

G(r,r',w) = G{r ~ r' ,uj) + G{r,u)f{u)G{~r' ,uj), (10) 

with the impurity induced T matrix can be obtained as'^?, 

^^"^ ^ l-G(Ur3 ^"--' (^^^ 



where G(w) = {l/N)Y,kG{k,uj). In this case, LDOS of 
the system can be obtained as, 

P{^,^) = --ImG(r,w) = pq{u) +6p{r,uj), (12) 

TT 

where the homogeneous density of states po(w) = 
— ImG(0, a;)/7r, and is uniform in real-space, therefore it 
reflects a homogenous background, while 6p{r,uj) is the 
modulation for the homogeneous density of states due to 
the presence of the impurity scattering potential ([9]) , and 
can be obtained as, 

(5p(r,w) = pir,u}) - poiuj) 

= --lm[G{r,Lj)f{Lj)G{-r,uj)]n, (13) 

TT 

and its Fourier transform is evaluated explicitly as, 
Sp{q,u) = -ilm (j^Y.^Gik + q,u)f{u)G{Kuj)li^ .(14) 
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This LDOS is closely related to the differential conduc- 
tance g{r,uj) since the result of the differential conduc- 
tance g{r,uj) is proportional to /9(r, oj) at location r and 
energy uj—eY. However, the intense atomic-scale spatial 
fluctuations in electronic structure cause systematic er- 
rors in setting the STM tip elevation for the experimental 
measurement of the differential conductance g(r, w)i2r— , 
therefore the ratio of differential conductances at oppo- 
site bias in Eq. ^ or its equivalent, the ratio of LDOS 
at opposite energy. 



(a) ra/J=-0.01 



(d) (i)/J=-0.01 



p(r,-w)' 



(15) 



has been measured experimentally for an enhancement 
of the QSI signaturesi^rii. Since the condition po{ijj) ^ 
Sp(r, uj) is well satisfied for cuprate superconductors, 
then the ratio of LDOS at opposite energy can be 
obtained approximately as^. 
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where Zq{uj) — pq{uj)/po{—lu), and only the first-order 
modulation Sp(r, ±uj) is kept, then the modulation of the 
ratio of LDOS at opposite energy can be expressed as. 



SZ{r,u) = Zir,Lu)-Zoiio) 
'Sp{r,uj) 



Zo{uj) 



Ml", -w) 



Po(w) 
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and its Fourier transformation is evaluated explicitly as. 



SZiq,Lj) ^ Zoiio) 



Spiq,(^) _ Sp{q, 
Po(w) 



poi-uj) 



(18) 



It has been arguedi^rJl that the observed QSI from 
SZ{r, uj) or 6Z{c{, lu) is an intrinsic phenomenon free from 
any set-point-rclated issues, which inevitably contami- 
nate g{r,uj) and ^(q, w). 

We are now ready to discuss the QSI phenomenon in 
cuprate superconductors in the SC state and its extinc- 
tion at the high energy. In cuprate superconductors, al- 
though the values of J and t is believed to vary somewhat 
from compound to compound, however, as a qualitative 
discussion, the commonly used parameters in this paper 
are chosen as t/J = 2.5 and t'/t — 0.3. In this case, we 
have performed a calculation for the Fourier transformed 
LDOS (5p(q, u) in Eq. (IT4l) and the Fourier transformed 
ratio of LDOS at opposite energy 6Z{q,uj) in Eq. (IT51) . 
and the results of |(5/9(q, w)| as a function of momentum 
in the full Brillouin zone with energy (a) uj = — O.OIJ, 
(b) UJ ~ — 0.03J, and (c) uj = — 0.05J at temperature 
T = 0.002J for the doping concentration 6 = 0.15 in 
the presence of single point-like potential scatterer of 
strength V ~ O.IJ are plotted in Fig. 1 (left panel). 
For comparison, the corresponding results of |(5Z(q, w)| 
are also plotted in Fig. Id-f (right panel). It is shown 
that the results obtained from SZ{q,uj) are almost the 
same as those of (5p(q, w), reflecting an experimental fact 
that the measurement data from the ratio of LDOS at op- 
posite energy retain all the main information of QSI ob- 
served from LDOS. Moreover, both results from SZ{q, u) 
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FIG. 1: (Color online) The Fourier transformed LDOS as a 
function of momentum in the full Brillouin zone with energy 
(a) UJ = -O.OIJ, (b) UJ = -0.03J, and (c) uj = -0.05J at tem- 
perature T = 0.002J for the doping concentration 5 = 0.15 in 
the presence of single point-like potential scatterer of strength 
V = O.IJ (left panel), (d) to (f) are the corresponding results 
of the Fourier transformed ratio of LDOS at opposite energy 
(right panel). 



and (5/9(q, uj) are clearly fourfold symmetric and display 
numerous local maxima (bright regions) at different q for 
different energies, where all of these q vectors are consis- 
tent with the prediction from the octet model^i^. These 
bright regions in momentum-space display the intensity 
of the modulation for LDOS (Fig. la-c) or the intensity 
of the modulation for the ratio of LDOS at opposite en- 
ergy (Fig. Id-f). In particular, the bright regions near 
the center of each figure in Fig. 1 a-c or Fig. Id-f reflect 
long- wavelength inhomogeneity in 6p{q,ui) or SZ{q,uj), 
and are induced by the weak impurity scattering poten- 
tial All these results are in qualitative agreement 
with the STM experimental results^i^. 

To analyze this peak feature in Fig. 1 more clearly, 
we have made a series of calculations for the momen- 
tum dependence of (5p(q, uj) and 6Z{q, uj) with different 
energies, and the results of \6p{q,uj)\ as a function of 
momentum along the (a) [0, 0] [tt, tt] (nodal) direction 
and (b) [0,0] — > [tt, 0] (antinodal) direction with energy 
UJ = -O.OIJ, UJ = -0.02J, UJ = -0.03J, UJ = -0.04J, 
and UJ ~ — 0.05J (from bottom to top) at temperature 
T — 0.002 J for the doping concentration 5 = 0.15 in 
the presence of single point-like potential scatterer of 
strength V = O.IJ are plotted in Fig. 2 (left panel). 
For comparison, the corresponding results of |(5Z(q, w)| 
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FIG. 2: The Fourier transformed LDOS as a function of 
momentum along tiie (a) [0, 0] — !> [tt, n] direction and (b) 
[0, 0] [tt, 0] direction with energy cj = -O.OIJ, u = -0.02 J, 
uj = — 0.03J, uj — — 0.04J, and uj — — 0.05J (from bottom to 
top) at temperature T = 0.002 J for the doping concentration 
5 = 0.15 in the presence of single point-like potential scatterer 
of strength V = O.IJ (left panel), (c) and (d) are the corre- 
sponding results of the Fourier transformed ratio of LDOS at 
opposite energy (right panel). 
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FIG. 3: The Fourier transformed ratio of LDOS at opposite 
energy as a function of momentum along the (a) [0, 0] — )■ 
[tt, tt] direction and (b) [0, 0] [tt, 0] direction with energy 
cj = — 0.02J at temperature T = 0.002J for the doping con- 
centration (5 = 0.08, 5 — 0.11, and 5 = 0.15 (from bottom to 
top) in the presence of single point-like potential scatterer of 
strength V = O.IJ. 



are also plotted in Fig. 2c-d (right panel). In correspond- 
ing to the results in Fig. 1, the results in Fig. 2 show 
that although the peak in 6p{q, oj) (then the intensity of 
the modulation for LDOS) or the peak in (5Z(q, oj) (then 
the intensity of the modulation for the ratio of LDOS at 
opposite energy) is changing rapidly and intricately with 
energy, the peaks, whose q-vectors oriented towards the 



antinodal direction, appear at finite |q| at very low en- 
ergy and then move steadily inwards toward the center, 
while the peaks with q-vectors along the nodal direction, 
appear and move steadily to large |q| with increasing en- 
ergy, in qualitative agreement with the STM experimen- 
tal results^. However, although the position of the peak 
obtained from (5p(q, cj) is almost the same as that from 
(5Z(q, w), the peak intensity obtained from (5Z(q, w) is 
much stronger than that appeared in 5p{q, lo), and there- 
fore there is an enhancement of the QSI signatures in 
SZ{q,uj). This expected difference of the peak intensity 
obtained from (5/o(q, w) and SZ{ci,uj) can be understood 
from the definition of SZ{q, uj) in Eq. p8|) . Although the 
modulations in Sp{r, cu) and Sp(r, —oj) are spatially quite 
similar, the spatial-phase relation for QSI is not known 
precisely. However, SZ(q,uj) is sensitive to the relation 
signs between i5p(q, w) and 6p{c{, —lo), namely, taking the 
ratio of LDOS reduces the in-plane component^- . This is 
why the the QSI extinction at the high energy in cuprate 
superconductors in the SC state is firstly observed from 
the experimental measurement of 52'(q, w)i2rii. 

For a better understanding of the physical proper- 
ties of the QSI phenomenon in cuprate superconduc- 
tors, we have further performed a calculation for SZ(q, uj) 
with different doping concentrations, and the results of 
\SZ{q_,uj)\ as a function of momentum along the (a) 
[0,0] ^ [7r,7r] (nodal) direction and (b) [0,0] [tt, 0] 
(antinodal) direction with energy uj = —0.02 J at temper- 
ature T = 0.002J for the doping concentration S — 0.08, 
6 = 0.11, and 6 = 0.15 (from bottom to top) in the pres- 
ence of single point-like potential scatterer of strength 
V = O.IJ are plotted in Fig. 3. Obviously, the peak is 
doping dependent. With increasing the doping concen- 
tration, the weight of the peak increases. Furthermore, 
the position of the peak along the nodal direction ap- 
pears at finite |q| in the low doping concentration, and 
then moves towards to the center of the Brillouin zone 
with increasing doping concentration. In contrast to the 
case along the nodal direction, the position of the peak 
along the antinodal direction is located at finite |q| in the 
low doping concentration, and then is shifted to large |q| 
with increasing doping concentration. 

The energy dependence of the peaks in Fig. 2 has been 
used to extract the dispersion of the peaks, and the re- 
sults of (a) the positions and (b) amplitudes of the lowest 
energy peaks in |(5Z(q, w)| as a function of energy along 
the [0,0] ^ [7i',7r] direction for temperature T = 0.002 J 
at the doping concentration S = 0.15 in the presence of 
single point-like potential scatterer of strength V = O.IJ 
are plotted in Fig. 4. For comparison, the corresponding 
STM experimental results^^ of cuprate superconductors 
are also shown in Fig. 4 (inset). Our results show clearly 
that the sharp peak persists in a large momentum space 
region along the nodal direction and are in qualitative 
agreement with the STM experimental result- How- 
ever, the amplitude of the peak increases with increasing 
energy at the low energy, and reaches a maximum at 
the intermediate energy, then decreases rapidly with fur- 
ther increasing energy. Apparently, there is a substantial 
difference between theory and experiment, namely, the 
amplitude of the peak obtained from theory does not di- 
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FIG. 4: (a) The positions and (b) amplitudes of the lowest 
energy peaks in the Fourier transformed ratio of LDOS at op- 
posite energy as a function of energy along the [0, 0] [n, n] 
direction for temperature T — 0.002J at the doping concen- 
tration 5 — 0.15 in the presence of single point-like potential 
scatterer of strength V = O.IJ. Insets: the corresponding 
experimental results for cuprate superconductors taken from 
Ref M 
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FIG. 5: The amplitudes of the lowest energy peaks in the 
Fourier transformed ratio of LDOS at opposite energy as a 
function of energy along the [0, 0] — >■ [vr, tt] direction with tem- 
perature T = 0.002 J at the doping concentration 5 = 0.15 for 
parameters F = 0.005J and a = 0.2 in the presence of single 
point-like potential scatterer of strength V = O.IJ. Inset: the 
corresponding experimental result for cuprate superconduc- 
tors taken from Ref. 



minish to zero at the high energy. This obvious weakness 
is a natural consequence of the neglect of the imaginary 
part of the self-energy (then the scattering rate). To 
obtain the full charge carrier Green function ([6]) within 
the kinetic energy driven SC mechanism, we^^'^'^ have 
made a self-consistent calculation in the static limit ap- 
proximation for the real part of the charge carrier self- 
energy Ke'E\'^\uj) induced by the interaction between the 
charge carriers and spins, where its antisymmetric part 

has been obtained as Zf^p = 1 — ReI]J^^(k, w = 0) |k=[;T,o] 
and therefore is closely related to the charge carrier 
quasiparticle coherent weight, while its symmetric part 
ReS^g-* (k, cj = 0) is a constant around the Fermi surface, 
and it just renormalizes the chemical potential. How- 
ever, the imaginary part of the charge carrier self-energy 
lml]['^\'k, uj) has been dropped^^'^^. This leads to a fact 
that the imaginary part of the renormalized electron self- 
energy in the dressed electron Green's function (ITUl) due 
to the presence of the impurity scattering potential © 



has been neglected. In this case, for the present discus- 
sions of the QSI extinction at the high energy, we need 
to introduce an imaginary part of the electron self-energy 
(then the scattering rate) in the electron Green function 
as, 



zlml](w) = i(T + auj), 



(19) 



for a compensation due to the neglect of the imaginary 
part of the renormalized self-energy, where F and a are 
constant. It has been argued^ that this type of the en- 
ergy dependent inelastic scattering rate ([TO]) seems to be 
an intrinsic property of the electronic structure of cuprate 
superconductors, since it has been used to provide a ro- 
bust fit for the spatially inhomogeneous differential con- 
ductances of cuprate superconductor based on the phe- 
nomenological d-wave BCS formalism. By considering 
this scattering rate in the electron Green function, 
we therefore find that the dispersing QSI always disap- 
pears at the high energy. To show this point clearly, 
we have made a series of calculations for SZ(q,u}) with 
different momenta and energies, and the result of the 
amplitudes of the lowest energy peaks in |(5Z(q, a;)| as 
a function of energy along the [0,0] [t^jT^] direction 
for temperature T = 0.002J at the doping concentra- 
tion S = 0.15 for parameters F = 0.005 J and a = 0.2 
in the presence of single point-like potential scatterer of 
strength V — O.IJ is plotted in Fig. 5 in comparison with 
the corresponding STM experimental result^ of cuprate 
superconductors (inset). In comparison with the result 
in Fig. 4b, although the amplitude of the peak is severely 
suppressed due to the presence of the scattering rate ([19]), 
it diminishes to zero (then the QSI extinction) at the high 
energy. Moreover, during the calculations, we find that 
the constant scattering rate F in Eq. (flQl) only plays 
a subsidiary role, while the QSI modulation is further 
weakened by the effective scattering rate acd. The similar 
conclusion has been obtained in Ref. [s^ where they also 
find that the effective scattering rate aui plays a key role. 
Using a reasonably estimative value of J ~ 120 ~ 150 
meV, the peak disappears around the energy 6 meV'^7.5 
meV. Apparently, there is still a substantial difference be- 
tween theory and experiment, namely, the energy value 
of the QSI extinction calculated theoretically is smaller 
than the corresponding energy value measured in the ex- 
periment. However, upon a closer examination one can 
see immediately that the main difference is due to fact 
that the calculated peak energy decreases rapidly with 
energy at high energy. This energy value of the QSI ex- 
tinction is also smaller than the corresponding value of 
the charge carrier pairing gap2^. The simple t-J model 
can not be regarded as a complete model for the quantita- 
tive comparison with cuprate superconductors, however, 
as for a qualitative discussion in this paper, the overall 
shape seen in the theoretical result is qualitatively con- 
sistent with that observed in the STM experiments^. 

A nature question is why QSI below the high energy 
in cuprate superconductors and its extinction at the high 
energy can be described qualitatively in the framework of 
the kinetic energy driven SC mechanism in the presence 
of a single point-like impurity. This may be understood 
from the kinetic energy driven SC mechanism itselfii. 
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As we have mentioned above, in the framework of ki- 
netic energy driven SC mechanism, the electron Cooper 
pairs originating from the d-wave charge carrier pairing 
state are due to the charge-spin recombination, therefore 
there is a coexistence of the electron Cooper pair and 
short-range AF correlation, and hence the short-range 
AF fluctuation can persist into the SC state^*". In par- 
ticular, this charge-spin recombination is characterized 
by a convolution of the spin Green's function and charge 
carrier Green's functionii. The main consequence of this 
convolution is that the hole-like charge carrier d-wave 
BCS formalism is transferred into the electron-like d- 
wave BCS formalism^^. In other words, main difference 
between the hole-like charge carrier Green's function ^ 
and the electron-like Green's function Q is a shift of 
the momentum by the AF wave vector Q = [tt, vr] in the 
Brillouin zone-^. This leads to a fact that the quasi- 
particle states near the antinodal points are broadened. 
Moreover, these quasiparticle states around the antin- 
odal regime are further suppressed due to the presence 
of the impurity scattering potential ([9]). On the other 
hand, the kinetic energy driven d-wave SC state is con- 
trolled by both the SC gap function and the SC quasipar- 
ticle coherent weight, which indicates that only coherent 
Bogoliubov quasiparticles at the Fermi surface are avail- 
able for superconductivity since everything happens at 
the Fermi surface. In this case, although the Bogoliubov 
quasiparticles can be excited by either positive or nega- 
tive energies^, the excitation spectrum is a particle-hole 
symmetric thin Dirac cone around the nodal regime, then 
QSI below the high energy and its extinction at the high 
energy occur in this case, where this high energy which 
sets the QSI extinction is determined not only by the 
d-wave SC gap function alone but also by the coherent 
Bogoliubov quasiparticles located on the thin Dirac cone. 

In summary, within the framework of the kinetic en- 
ergy driven SC mechanism, we have discussed the QSI 
phenomenon in cuprate superconductors in the SC state 
by considering the single point-like impurity scattering 
potential. This QSI is characterized by the peaks in 
LDOS or the ratio of LDOS at opposite energy. By 
calculation of the momentum and energy dependence of 
the Fourier transformed LDOS i5p(q, w) and the Fourier 
transformed ratio of LDOS at opposite energy 5Z{q,uj), 
we have shown that the remarkable QSI phenomenon ob- 
served from STM experiments on cuprate superconduc- 
tors can be described qualitatively by a single point-like 
impurity in the kinetic energy driven homogeneous d- 
wave SC state. The amplitude of the peak increases with 
increasing energy at the low energy, and reaches a maxi- 
mum at the intermediate energy, then diminishes to zero 
at the high energy. The theory also predicts that with 
increasing doping concentration, the position of the peak 
along the nodal direction moves towards to the center of 
the Brillouin zone, while the position of the peak along 
the antinodal direction is shifted to large momentum re- 
gion, which should be verified by further STM experi- 
ments. 

Within the framework of the kinetic energy driven SC 
mechanismii, our present results of QSI below the high 



energy and its extinction at the high energy due to the 
presence of the impurity scattering potential © are qual- 
itatively consistent with the recent STM experimental 
data^ii^. Establishing this agreement is important to 
confirming the nature of a single impurity in the kinetic 
energy driven homogeneous d-wave SC state in cuprate 
superconductors. Although a quantitative understanding 
of the QSI phenomenon is not straightforward to obtain, 
and depends rather sensitively on the nature of the scat- 
tering medium, in this paper we are primarily interested 
in exploring the general notion of the QSI phenomenon 
induced by a single impurity in the kinetic energy driven 
homogeneous d-wave SC state. The qualitative agree- 
ment between the present theoretical results and STM 
experimental data also show that the presence of im- 
purities plays a crucial role for the QSI phenomenon in 
cuprate superconductors in the SC state. Although this 
QSI extinction at the high energy^'* also can be fitted by 
using a phenomenological d-wave BCS formalism with 
the same imaginary part of the self-energy (jl9p . how- 
ever, no explicit calculations of the doping dependence 
of |i5Z(q,a;)|, as shown in Fig. 3, can be made within a 
phenomenological d-wave BCS formalism. 

Under the kinetic energy driven SC mechanism, the 
external magnetic field aligns the spins of the unpaired 
electrons, then the singlet charge carrier pairs can not 
take advantage of the lower energy offered by a spin- 
polarized state. In this case, the magnetic field depen- 
dence of the penetration depth and superfluid density 
have been studied recently'^®, where the superfluid den- 
sity decreases with increasing magnetic field, in agree- 
ment with the experimental results'^^. With this study"^-, 
we therefore can predict that QSI should be magnetic 
field dependent, and then the energy value of the QSI 
extinction decreases with increasing magnetic field. One 
of the typical features of d-wave superconductivity in 
cuprate superconductors is the particle-hole symmetric 
octet of dispersive Bogoliubov quasiparticle interference 
modulations^"—. However, the recent STM experiments^ 
have reported on this octet's evolution from low temper- 
atures to well into the underdoped pseudogap regime, 
where no pronounced changes occur in the octet phe- 
nomenology at the superconductor's critical temperature 
Tc, and it survives up to at least temperature T 1.5Tc. 
Thus an important issue is how this octet's evolution 
from the SC state into the pseudogap regime is fitted 
within the framework of the kinetic energy driven su- 
perconductivity. These and the related issues are under 
investigation now. 
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